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"Yale University, New Haven, Connecticut 06511, USA 
(Dated: February 7, 2008) 

We report the first measurement of the branching fraction /oo for Y(4S) — » B°B° . The data 
sample consists of 81.7 fb _1 collected at the T(4S) resonance with the BABAR detector at the PEP- 
II asymmetric-energy e + e~ storage ring. Using partial reconstruction of the decay B° — > D* + l~ui 
in which only the charged lepton and the soft pion from the decay D* + — » D°tv + axe reconstructed, 
we obtain /oo = 0.487 ± O.OlO(stat) ± 0.008(sys). Our result does not depend on the branching 
fractions of B° — > D* + l~T>i and D* + — > D°tt + decays, on the ratio of the charged and neutral B 
meson lifetimes, nor on the assumption of isospin symmetry. 

PACS numbers: 14.40.Nd, 13.20.He 



Isospin violation in the decay T(4S) — » BB will lead 
to a difference between the branching fractions /oo = 
B{T{AS) -► B°B°) and /+_ ee B{T{AS) -* B + B~). Pre- 
dictions for the ratio = /_| If on range from 1.03 to 

1.25 Q. Measurements of R+l° lllllllSll have been 
made assuming isospin symmetry in specific decay rates 
and resulting in an average value of 1.006±0.039 con- 
sistent with isospin conservation in T(4S I ) decays to BB. 
To date no measurement has been made of either /oo or 

/-I In this paper we report the first direct measurement 

of /oo • It is completely independent of the previous mea- 
surements of Independent measurements of /oo 
and can be used to constrain the T(AS) — ► non- BB 
fraction. The /oo value is important for measuring ab- 
solute T(AS) branching fractions and for measuring V c b, 
the Cabibbo-Kobayashi-Maskawa matrix element. 

The data sample used in this analysis consists of 
81.7fb _1 collected at the T(45) resonance (on-resonance) 
and 9.6 fb" 1 collected 40MeV below the resonance (off- 
resonance). The on-resonance data sample has a mean 
energy of 10.580 GeV and an energy rms spread of 
4.6 MeV. Due to the small spread, any plausible energy 
dependence of /oo has a negligible effect on the central 
value. A simulated sample of BB with integrated lumi- 
nosity equivalent to approximately three times the data 
is used for background studies. 

A detailed description of the BABAR detector and the 
algorithms used for track reconstruction and particle 
identification is provided elsewhere Q. A brief sum- 
mary is given here. High-momentum particles are recon- 
structed by matching hits in the silicon vertex tracker 
(SVT) with track elements in the drift chamber (DCH). 
Lower momentum tracks, which do not leave signals on 
many wires in the DCH due to the bending induced by 
a magnetic field, are reconstructed in the SVT alone. 
Electrons are identified by the ratio of the track momen- 
tum to the associated energy deposited in the calorime- 
ter (EMC), the transverse profile of the shower, the en- 
ergy loss in the drift chamber, and information from 
a Cherenkov detector (DIRC). Muons are identified in 
the instrumented flux return (IFR), composed of resis- 
tive plate chambers and layers of iron. Muon candidates 
are required to have a path length and hit distribution 
in the IFR and energy deposition in the EMC consis- 



tent with that expected for a minimum-ionizing particle. 
The BABAR detector Monte Carlo simulation is based on 
GEANT4 0. 

We select the decays B° -> D*+£-D e , D*+ -► D°ir+ 
(£ = e, fi). The inclusion of charge-conjugate reactions 
is implied throughout this paper. The sample of events 
in which at least one B° — > D* + £~i?i candidate decay is 
found is labeled the "single-tag sample" . The number of 
signal decays in this sample is 

N s = 2N BB -f m e a B(B° -> D*+rP e ), (1) 

where N BB is the total number of BB events in the 
data sample and e s is the reconstruction efficiency for 
B° — ► D* + £~D ( . We determine N BB = 88.7 million events 
by counting the number of hadronic decays in the on- 
resonance data and subtracting the e + e~ — * qq (q = u, 
d, s, or c quark) component using off-resonance data, 
as described in detail in Ref. 0. The error in N BB is 
1.1% and is dominated by systematic uncertainties. We 
attribute all BB pairs to T(AS) decays. 

The number of signal events in the subset in which 
two B° — > D* + £~i>i candidates are found is labeled the 
"double-tag sample" . The number of such events is 

N d = N BB f 0Q e d [B(B° -> D*+£-9e)] 2 , (2) 

where e c i is the efficiency to reconstruct two B° — > 
D* + £~Di decays in the same event. From Eq. (JJJ and 
Eq. ®, /oo is given by 

CN 2 

/0 ° = 4AVV^' (3) 

where we have defined C = Sd/e 2 . The value of C is 1 if 
the efficiencies for detecting each B meson are uncorre- 
cted in double-tag events, which, given the pseudoscalar 
nature of B mesons and the proximity of the Y(AS) to 
the BB threshold, is expected. Using the Monte Carlo 
simulation we determine C = 0.995 ± 0.008, where the 
error is due to the finite size of the simulated sample. 

We select the decays B° — > D* + £~vi with a partial 
reconstruction technique 0,0, 1X3- E3| - In this technique, 
only the lepton from the decay B° — > D* + £~vt and the 
soft pion from the decay D* + — » D°ir + are reconstructed. 
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No attempt is made to reconstruct the D°, resulting in 
a high reconstruction efficiency. 

The B° decay point is determined from a vertex fit 
of the soft-pion and lepton tracks, with the vertex con- 
strained to the beam spot position in the x — y plane. We 
only use events with vertex-fit probability, Vy, greater 
than 0.1% to optimize a signal-to-background ratio. 

We select hadronic events by requiring at least four 
charged particle tracks reconstructed in the SVT and the 
DCH. To reduce non-BB background, the ratio of the 
second to the zeroth Fox- Wolfram moments 14], i?2 = 
H2/H0, is required to be less than 0.5. 

To suppress leptons from charm decays, all lepton can- 
didates are required to have momenta between 1.5 GeV/c 
and 2.5GeV/c in the e + e~ center-of-mass frame. Soft 
pion candidates are required to have center-of-mass mo- 
menta between 60MeV/c and 200MeV/c. As a conse- 
quence of the limited phase space available in the D* + 
decay, the soft pion is emitted nearly at rest in the D* + 
rest frame. The D* + four-momentum can therefore be 
computed by approximating its direction as that of the 
soft pion, and parameterizing its momentum as a linear 
function of the soft-pion momentum, with parameters 
obtained from a Monte Carlo simulation. The presence 
of an undetected neutrino is inferred from conservation 
of momentum and energy. The neutrino invariant mass 
squared is calculated as 



M 2 ee (E h 



cam 



E 



D> 



Et) 2 - (pd* + Pi 



(4) 



where -E^eam ^ s na ^ ^he center-of-mass energy and 
E e (E D *) and p e (pz?*) are the center-of-mass energy 
and momentum of the lepton (the D* meson). We set 
Pb = 0, which introduces a negligible spread in M 2 
compared with the approximation of the D* momen- 
tum based on the soft pion. For signal decays that are 
properly reconstructed, the Ai 2 distribution peaks near 
zero. Background events, however, are spread over a 
wide range of M 2 values. We define a signal region 
(M 2 > -2 GeV 2 /c 4 ) and a sideband region (-8 < M 2 < 
-4 GeV 2 /c 4 ). 

We use the symbol AI 2 to denote M 2 for any candi- 
date in the single-tag sample. In the double-tag sam- 
ple, we randomly choose one of the two reconstructed 
B° — > D* + £~Dg candidates as "first" and the other as 
"second". Their M 2 values are labeled M\ and M\, re- 
spectively. We require that Ai 2 fall in the signal region. 

The single-tag and double-tag samples have several 
types of background: continuum, combinatorial BB, and 
peaking BB. The combinatorial BB background orig- 
inates from random combinations of reconstructed lep- 
tons and soft pions. The peaking BB background is com- 
posed of B — » D*tt£D£ decays with or without an excited 
charmed resonance D** [J^, where the reconstructed soft 
pion comes from the decay D* + — » D°tt + , leading to an 
accumulation of these events at high values of Ai 2 . The 
peaking BB background is suppressed by the requirement 



pe > 1.5GeV/c on the lepton center-of-mass momentum. 
Such events have an Ai 2 distribution that is different 
from the signal, allowing us to extract their contribution 
in the signal region. 

The double-tag sample contains two additional types of 
background: events in which the first candidate is com- 
binatorial background and the second is signal (called 
Ai ^-combinatorial background) and events in which the 
first candidate is peaking background and the second is 
signal (called Al^-peaking background). 

To determine N s and Nd, we perform binned x 2 fits 
to one-dimensional histograms of the Ai 2 and Ai 2 dis- 
tributions of on-resonance data events, ranging from —8 
to 2 GeV 2 /c 4 . Before fitting, we subtract the continuum 
background contribution from the histograms. This is 
determined using the AI 2 and Ai 2 , distributions of off- 
resonance data, scaled to account for the ratio of on- 
resonance to off-resonance luminosities and the center- 
of-mass energy dependence of the continuum production 
cross-section. In addition, the contributions of the Ai\- 
combinatorial (3%) and A( 2 -peaking (1%) backgrounds 
are subtracted from the Ai\ histogram before doing the 
fit. The contribution of the A4 2 -combinatorial back- 
ground is determined from sideband data. The Ai\- 
peaking background is determined with simulated events. 

After the subtraction, the Ai 2 and Ai 2 histograms are 
fit separately, to a function whose value for bin j of the 
histogram is 



(5) 



where A* is the number of events of type t (t = signal, 
combinatorial, peaking) populating the histogram, and 
Pj is the bin j value of a discrete probability density 
function (PDF) obtained from simulated events of type 
t, normalized such that J2j -Pj — !• The fit determines 
the parameters A* by minimizing 



E 

3 



(6) 



'Hi 



where Hj is the number of entries in bin j of the data 
histogram being fit; erg. is the statistical error on Hj, in- 
cluding uncertainties due to the background subtractions 
described above; and af j is the error on fj, determined 
from the errors on Pj , which are due to the finite size of 
the simulated sample. 

The results of the fits are presented in Table|U The Ai 2 
and AI 2 distributions are shown in Fig. GJ The fits yield 
N s = 786200 ± 1900 (Confidence Level (C.L.)= 11%) 
and N d = 3560 ± 70 (C.L. = 82%). Equation then 
gives foo = 0.487 ± 0.010, where the error is due to data 
statistics. 

To determine how well the simulation reproduces the 
AI 2 and A4 2 distributions for the combinatorial back- 
ground in the data, we study the distributions for a 
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TABLE I: Numbers of entries of different types in the M 2 S 
and _M| histograms in the signal region. 



Source 


Ml 


Ml 


Signal 


786200 ± 1900 


3560 ± 70 


Combinatorial BB 

Peaking BB 

Continuum 

M i-combinatorial 

A^f-peaking 


558080 ± 470 
68170 ± 260 
240600 ± 1400 


1510 ±20 
300 ± 20 
160 ± 40 
180 ± 20 
60 ± 10 



sample of same-charge candidates, in which the lepton 
and soft pion have the same electric charge. We fit the 
continuum-subtracted M? 3 and Mil histograms of the 
same-charge sample using the function /j = NP^, where 
Pj is the bin j value of the PDF of same-charge simulated 

BB events, normalized such that J2j P'j = 1' an d the pa- 
rameter N is determined by the fit. The histograms, 
overlaid with the fit function, are shown in Fig. [5] The 
accumulated differences D = J2j(Hj ~ f'j) between the 
same-charge data histograms H'^ and the fit functions are 
summarized in Table [n] Their consistency with zero in- 
dicates that the distributions of simulated combinatorial 
BB background events do not lead to significant fake sig- 
nal yields. Nevertheless, we evaluate a systematic uncer- 
tainty on the modeling of the combinatorial background 
based on the observed difference in the like-sign sample. 



TABLE II: The difference D = J2j( H 'j ~ f'j) between the 
same-charge data histogram and the fit function, summed 
over the signal region or over the whole region of the Ml 
and M\ distributions. 



Fit 


Signal re; 


;ion 


Whole region 


parameter 


M 2 


Ml 


Ml Ml 


D 


-1300 ±2100 - 


-80 ±80 


700 ± 3000 70 ± 80 


C.L.(%) 


57 


78 


94 98 



We evaluate the absolute systematic uncertainties in 
/oo due to the M.\ -combinatorial subtraction (0.0005), 
the .Mf-peaking background (0.0005), the value of C 
due to the track multiplicity dependence of the efficiency 
(0.0015), the finite size of the simulated sample (0.002), 
the same-charge sample (0.0025), the impact of a pos- 
sible contribution of non- BB decays of the T(45) 03 
(0.0025), the peaking background composition (0.004), 
and the total number of BB, N bS (0.0055). 

The dominant contribution to the systematic error 
comes from a 1.1% systematic uncertainty in N BB -, due 
mainly to the uncertainty in the tracking efficiency. The 
peaking BB background is estimated from the simulated 
sample containing all D** resonances and non-resonant 
events. We vary the ratio of the branching fraction of the 
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FIG. 1: The Ml (top) and M% (bottom) distributions for 
the on-resonance sample. The continuum background has 
been subtracted from the distributions. In addition, the M\- 
combinatorial and the A^i-peaking backgrounds have been 
subtracted from the M\ distribution. The levels of the simu- 
lated signal, peaking BB, and combinatorial BB background 
contributions are obtained from the fit. 



resonant and the non-resonant production such that the 
variation of this ratio is wide enough to include poorly 
known decays. We repeat the analysis procedure to de- 
termine N s and Nd- The uncertainties due to the lep- 
ton and soft-pion momentum spectra are negligible. We 
combine the uncertainties given above in quadrature to 
determine an absolute systematic error of 0.008 for /oo- 
In summary, we use a partial reconstruction of the de- 
cay B° — > D* + t~Di to obtain the result 



/oo = 0.487 ± O.OlO(stot) ± 0.008(sys), 



(7) 
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FIG. 2: The M 2 S (top) and Ml (bottom) distributions for 
the same-charge on-resonance sample. The continuum back- 
ground has been subtracted from the distributions. The M\- 
combinatorial and the A'ff-peaking backgrounds have been 
subtracted from the M\ distribution. The level of the simu- 
lated combinatorial BB background is obtained from the fit. 



where the first error is statistical and the second is sys- 
tematic. This result is the first, precise, and direct mea- 
surement of foQ. Since this measurement is made by com- 
paring the numbers of events with one and two recon- 
structed B° — » D* + £~&e decays, it does not depend on 
branching fractions of B° — > D* + l~Vi and D* + — > D°ir + 
decays, on the ratio of the charged and neutral B meson 
lifetimes, nor on the assumption of isospin symmetry. By 
combining our /oo result with the world average of 
noted in the introduction, we add the errors quadratically 
to obtain /+_ = 0.490 ±0.023. Thus we find the fraction 
of Y(AS) -> non-55 to be 1 - /oo - /+- = 0.023 ± 0.032. 
If /oo + /h — = 1, our /oo result can be averaged with 
R+/° to yield / 00 = 0.494±0.008, /+_ = 0.506±0.008, 
and /+_//oo = 1-023 ± 0.032. This value of /+_//oo is 
in good agreement with isospin conservation in T(45*) — > 
BB within errors. 
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